Multi-epitope vaccine design using in silico analysis of glycoprotein and nucleocapsid of NIPAH virus

According to the 2018 WHO R&D Blueprint, Nipah virus (NiV) is a priority disease, and the development of a vaccine against NiV is strongly encouraged. According to criteria used to categorize zoonotic diseases, NiV is a stage III disease that can spread to people and cause unpredictable outbreaks. Since 2001, the NiV virus has caused annual outbreaks in Bangladesh, while in India it has caused occasional outbreaks. According to estimates, the mortality rate for infected individuals ranges from 70 to 91%. Using immunoinformatic approaches to anticipate the epitopes of the MHC-I, MHC-II, and B-cells, they were predicted using the NiV glycoprotein and nucleocapsid protein. The selected epitopes were used to develop a multi-epitope vaccine construct connected with linkers and adjuvants in order to improve immune responses to the vaccine construct. The 3D structure of the engineered vaccine was anticipated, optimized, and confirmed using a variety of computer simulation techniques so that its stability could be assessed. According to the immunological simulation tests, it was found that the vaccination elicits a targeted immune response against the NiV. Docking with TLR-3, 7, and 8 revealed that vaccine candidates had high binding affinities and low binding energies. Finally, molecular dynamic analysis confirms the stability of the new vaccine. Codon optimization and in silico cloning showed that the proposed vaccine was expressed to a high degree in Escherichia coli. The study will help in identifying a potential epitope for a vaccine candidate against NiV. The developed multi-epitope vaccine construct has a lot of potential, but they still need to be verified by in vitro & in vivo studies.


Introduction
Emerging/remerging infectious viral diseases that have posed a hazard to humanity throughout history [1].The ongoing COVID-19 outbreak has had a catastrophic impact on both the global economy and health [2].Numerous additional viruses, in addition to COVID-19, share this mode of transmission [3,4].In view of this pandemic, the international community must develop countermeasures in advance of virus outbreaks in order to be ready to combat them.
The Hendra virus (HeV) and Nipah virus (NiV), which had previously been responsible for severe and frequently fatal epidemics in both humans and animals, were discovered in Australia and Malaysia, respectively, in the 1990s.HeV is only present in Australia, but recent NiV outbreaks were initially discovered in Bangladesh, India, and the Philippines.According to several studies the mortality rate for the highly fatal zoonotic virus known as the NiV ranges from 40 to 91% [5][6][7][8].Based on the WHO R&D Blueprint list of priority diseases, it is recommended that the development of the NiV vaccine be expedited.According to stage III zoonotic disease classification standards [9,10], NiV infection can transmit to people and result in isolated epidemics.
The first-ever outbreak in southern India occurred in May 2018 in the state of Kerala, with a total of 19 NiV cases, 17 of which were fatal.Pteropus giganteus bats from the outbrake regions of the Kozhikode, Kerala, showed positivity for NiV.Since some NiV strains may spread from person to person and have a high mutation rate, there is an increased chance that NiV will spread globally and cause a pandemic [11].There is presently no vaccine or licensed therapy for the treatment of NiV.While these outbreaks often tend to be modest, we cannot rule out the possibility that the virus will evolve and/or spread.In addition to being severe, these outbreaks also have a negative influence on society.Infection could spread quickly and humanity could experience its most devastating pandemic if a NiV strain evolved to be human-adapted and infected communities in Southeast Asia, where pigs are a major export product and there are high densities of both humans and pigs [9,12,13].
NiV's genome is a non-segmented single-stranded RNA that is around 18.2 kilobases in size and encodes six structural proteins: nucleocapsid (N), glycoprotein (G), matrix protein (M), fusion protein (F), phosphoprotein (P), and RNA polymerase (L) [14].The NiV glycoprotein (NiV-G), is a type II transmembrane glycoprotein that is essential for the virus to enter host cells.A homotrimeric protein, the NiV-G protein is composed of three identical subunits, each of which is made up of two domains: an N-terminal receptor binding domain (RBD) and a C-terminal fusion domain (FD) [15].The RBD of the NiV-G virus is in charge of identifying and attaching to host cell receptors, which is what permits the virus to enter the cell.The cellular receptor for the NiV has been discovered to be the membrane-bound proteins ephrin-B2 and ephrin-B3, which play a crucial role in cell signalling and development.A recent investigation allowed for this discovery.The attachment of the NiV-G RBD to these receptors is made possible by its high affinity and specificity, and it is a vital stage in the process by which the virus enters host cells [16].The C-terminal fusion domain (FD) of NiV-G aids in the fusion of the host cell's membrane and viral membrane, allowing the viral entry in cytoplasm of the host cell's.The fusion process is initiated by conformational alterations in the NiV-G protein brought on by receptor contact [17].The very conserved NiV nucleocapsid protein, also known as NiV-NP, is essential for both viral assembly and replication.The multifunctional NiV-NP participates in viral transcription, RNA synthesis, and genome packing [18].The ribonucleoprotein complex that the NiV-NP forms serves as a template for viral transcription and replication.By attaching to the virus's RNA genome, this complex is created.In addition, the NiV-NP interacts with the L and P proteins of the viral polymerase complex to promote the synthesis of viral RNA.Additionally, as the virus is being put together, the NiV-NP is in charge of attracting the viral RNA genome into the developing viral particle.According to Wang, et al., this is one of the NiV-NP numerous significant functions in genome packing.It has been proven that antibodies that are focused on NiV-NP can successfully shield animals from NiV infection [19].According to Weingartl H and colleagues, the immune system also targets the NiV-NP [20].Immunoinformatic approaches are a powerful method for screening entire pathogen proteomes for potential immunogenic areas in this important race to create a vaccine [21].For the current investigation, we used an in silico methods to create a multiepitope vaccine against NiV.

Methodology
To identify potential vaccine candidates against NiV, a computational approach was used in the current investigation.In Fig 1 we see a flowchart that summarizes the methodology for using a reverse vaccinology technique to build an epitope-based vaccine.

Retrieval of NiV protein sequence and the analysis of the physiochemical and antigenic properties
From the NCBI database, the attachment glycoprotein (G) and nucleocapid (N) proteins of NiV were obtained.

Cytotoxic T lymphocytes (CTL) and helper T lymphocytes (HTL) epitope prediction
The Immune Epitope Database (IEDB), a comprehensive resource for epitope prediction that includes multiple algorithms and tools, including artificial neural networks and support vector machines.MHC-I and MHC-II epitopes for NiV-G and NiV-NP were predicted using the server's default parameters.Accordingly, the percentile score (lower the percentile value, higher the binding affinity of the predicted epitope) was chosen as <0.5 for both MHC-I and MHC-II alleles.

Immunogenicity prediction of MHC class I epitopes
In order to further analyse the immunogenicity of the MHC-I epitopes identified from the IEDB server, we used the IEDB MHC-I Immunogenicity tool, which examines and verifies exclusively 9-mer peptides [22].

Analysis of population coverage
The MHC HLA allele might vary depending on where in the world you are located due to the many different geographical regions.Because of this, conducting an analysis of population coverage is now required before building an effective vaccine construct.To evaluate the global coverage of allele-interaction epitopes of MHC-I and MHC-II, the IEDB population coverage analysis tool was used to conduct a population coverage study [24].The vast majority of possible prominent epitopes were found in various locations around the world.

Designing and physiochemical analysis of vaccine constructs
In silico vaccine constructs based on MHC-I, MHC-II, and B cell epitopes from NiV were created by combining the MHC-I, MHC-II, and B cell epitopes with AAY, GPGPG, and KK linkers, respectively [25,26].An EAAAK linker along with a cholera toxin adjuvant (Hou, J. et al. 2014; Kim, H. J., et al 2007) is also expected at the N terminus The protein sequence's physicochemical properties, such as molecular weight, aliphatic score, in-vitro half-life instability index, and grand average of hydropathy (GRAVY), were analyzed using the ProtParam program [27].

Prediction of conformational BCEs
The ellipro tool of IEDB was used to predict the conformational BCEs in the vaccine construct.Protein structures were uploaded in PDB format, while all other parameters were left at their default settings.Ellipro works by utilizing the geometric aspects of protein structure, including nearby cluster residues, the overall shape of the protein, and the remaining protrusion index [28].

Allergenicity and antigenicity
Based on the resemblance of potential vaccine candidates to known epitopes, AllerTOP v. 2.0 was used to evaluate the allergenicity of vaccine candidates [29].VaxiJen 2.0 server was used to evaluate the vaccine candidate's antigenicity [30].

Prediction of secondary structure
Psipred is an internet server that uses the amino acid sequences in a specific manner to produce the secondary structure of the multi-epitope construct (http://bioinf.cs.ucl.ac.uk/psipred/).This is one of the most popular servers for predicting protein secondary structures; it analyzes the PSI-BLAST output using two feed-forward neural networks [31].

Prediction, validation and refinement of the tertiary structure
I-TASSER, a server using the sequence-to-structure-to-function paradigm was used to predict the construct's tertiary structure [25].The created construct's 3D structure needs to be validated in order to look for any potential mistakes.The Ramachandran plot generated by PDBsum, a visual database that generated the tertiary structure of the vaccine construct, was used to validate it.

B-cell epitope (BCEs) prediction
The primary antigenic determinants recognized by the host immune system are BCEs, which describe a particular antigenic region that binds with the B lymphocytes.A vaccine candidate needs to start an efficient B-cell humoral immune response.By applying default parameters and a score of more than 0.8 for each epitope, we performed the artificial neural networkbased ABCpred server analysis uses recurrent neural network for the prediction of B cell epitopes and server has 65.93% accuracy [23].

Disulfide engineering of the vaccine construct
Disulfide engineering is a novel approach to enhance protein structure stability by introducing new disulfide bonds in the target protein through cysteine modification of protein structure residues.Hence, the Disulfide by Design 2 (DbD2) server was employed to detect residue pairings with mutational potential suitable for disulfide engineering [32].

Molecular docking & molecular dynamic simulation
Virus infection triggers an innate immune response through TLRs; in particular, it has been demonstrated that TLR can recognize the viral envelope protein.As a result, the target proteins against which vaccination candidates were developed were TLR3, TLR7, and TLR8.When predicting the binding energy between TLRs and vaccine constructs, the ClusPro 2.0 server was utilized [33].The experiment was carried out by uploading PDB files of the receptors and ligands into the server and submitting them with the default parameters.The best-scoring vaccine assembly complexed with TLRs was exposed to molecular dynamic simulation by iMODS, a online server after executing the protein-protein molecular docking [34], using the default parameters.Numerous antigen-presenting cells (APCs) including tissue dendritic cells and monocytes include the TLR-3 protein, which can activate the body's cells' unique recognition response to RNA virus infection.Predictions are made about protein complex dynamics simulation in terms of atomic B-factors, eigenvalue variance, deformability, elastic network, and covariance map.The capacity of each of a protein's residues to deform is largely responsible for a protein's deformability.The eigenvalue is associated to the energy needed to distort the specified structure; the lower the eigenvalue value, the easier the complex is to deform.Additionally, according to Sarkar et al, the eigenvalue of the protein complex in question determines the stiffness of its motion [35].

Analysis of the construct's immunological characteristics
The C-immSim immune server was used to model how an actual immune system would react to our final vaccine construct.Immune interactions and immunogenic epitopes can be simulated and predicted using this freely available web-server [36], which is based on the Position-Specific Scoring Framework (PSSM).The program was executed with the default settings, including a vaccination injection without lipopolysaccharide and time step of injection of 1, 42, and 84.

In silico cloning optimization of designed vaccine candidate
The multi-epitope vaccine candidate's codons were optimized, and reverse translations were examined using JCAT.The degree of expression of protein was calculated using a codon adaptation index (CAI).Additionally, JCat output included the proportion of GC content [37].The improved sequence was then cloned using DNASTAR into the pcDNA™3.1/V5-His-TOPO1expression vector.

Results
For the immunoinformatic investigation, NCBI was contacted for the attachment glycoprotein (G) and nucleocapid (N) proteins of NiV with accession numbers NP 112027.1 and NP 112021.1,respectively.

Cytotoxic T lymphocytes (CTL) and helper T lymphocyte (HTL) epitope prediction
The IEDB server was used to determine particular immunogenic CTL (9-mer) and HTL (15-mer) epitopes using the G and N protein sequences as input.The immunogenicity and percentile rank were used to filter and choose the predicted epitopes, and they were found to be likely antigens.For the various MHC I alleles of the G and N proteins, a total of 24 (Table 1) and 30 (Table 2) predicted epitopes were shown to be quite immunogenic and antigenic.The several CTL epitopes that were predicted exhibited strong binding affinity to multiple MHC I alleles.Thirteen HTL epitopes were G protein epitopes (Table 3), and 28 HTL epitopes for N protein, all predicted epiotpes were found to be IFN-inducing epitopes (Table 4).

Population coverage
It was anticipated that the vaccine that would be created would be able to protect a significant number of people all over the world.The evaluation of the distribution of CTL and HTL epitope alleles that were employed in the development of the vaccine construct was carried out with the assistance of the IEDB population coverage tool, which can be found at http://tools.percentages were as follows: 98.86% in Europe, 97.13% in the West Indies, 12.45% in Central America, 97.19% in North America, and 75.59% in South America.This proportion was 93.9 percent in Oceania.These findings show that the predictable epitopes are appropriate for use in vaccine construction and cover a wide spectrum of the human population.

Identification and prediction of B cell epitopes (BCEs)
The ABCpred server was utilized in order to make the prediction of the B cell epitopes using a cutoff value of 0.8.There were a total of 26 and 24 B cell epitopes predicted for the Nipah virus's G (Table 5) and N (Table 6) proteins, respectively.

Preparation of vaccine construct and physiochemical properties
Using the AAY, GPGPG, and KK linkers, the most promising MHC-I, MHC-II, and B cell epitopes were combined to create the vaccine design.The 124 amino acid long Cholera Toxin B (CTB), a potent adjuvant was used to enhance the immune response to co-administered antigens.As an adjuvant, CTB has been found to increase the production of antigen-specific antibodies and stimulate the activation of T-cells, leading to a more robust and effective immune

Secondary structure prediction
An accurate prediction of the secondary structure of the NiV vaccine was achieved using PSIPRED.Fig 4 shows the residues that make up the alpha helix to be represented in a pink color, the residues that make up the beta strand to be represented in a yellow color, and the residues that make up the coil to be represented in a gray color.Alpha helix is represented by the pink color, β strand by the yellow color, and coil by the gray color

Allergenicity and antigenicity
It was hypothesized that the multiepitope vaccine design would not trigger allergic reactions when run through the AllerTOPv2.0 algorithm.Our multiepitope construct received a predicted antigenicity score of 0.5512 from Vaxijen v2.0.This score was above the threshold value of 0.4%, suggesting that the construct possessed antigenic characteristics.

Conformational BCEs & disulfide engineering of the vaccine construct
A total of 5 conformational epitope has been predicted in the vaccine construct and have the higher score of 0.5 (Table 7).The vaccine construct's model identifies 32 residue pairs that may potentially create a disulfide bond, as predicted by the Disulfide by Design 2.13 server (Table 8).As per study by Craig DB, et al., for the formation of disulfide bonds the χ3 angle lies between -87˚and +97˚and energy value less than 2.2 kcal/mol, on the basis of the this only two amino acid pair (335-GLY-373-GLY & 24-GLN-28-ASP) were form the disulfide bond in the designed construct (Table 8).

Molecular docking of vaccine construct with receptors
Through molecular docking, the stability and binding affinity of a docked complex between a ligand and receptor molecule can be examined.For pathogen identification and immune response, toll-like receptor is a crucial human protein.In light of this, we select TLR3, 7, and 8 as the immunological receptors for the molecular docking.The revised 3D model of our finished vaccination is molecularly docked with the TLR3, 7, and 8 immunological receptors using the ClusPro 2.0 server (Kozakov D et al., 2017).The multi-epitope vaccination and TLRs (3, 7 and 8) were molecularly docked using the ClusPro 2.0 server, which created docked complexes with various cluster members and the lowest energy.All docked complexes were compared, and the Clusters with the lowest energy scores were chosen.The residues involved in the interaction between TLR3 (

Simulation of the vaccine-receptor complex using molecular dynamics
With the aid of the iMOD server, we run a molecular dynamics simulation to assess the stability and physical motions of the vaccine-TLR3, 7, and 8 docked complex (Lo ´pez-Blanco, J. R., Table 7. List of discontinuous epitopes and their predicted scores.Aliaga, 2014).Figs 7-9 illustrate the main-chain deformability for TLR 3, TLR 7, and TLR 8, respectively.Normal mode analysis indicates regions with high deformability when hinges are present, and B-factor values are proportional to root mean square values.Low RMSD and highly correlated areas in heat maps demonstrated the improved interaction of each residue.These findings demonstrate the structure of protein with its MNA mobility, the low deformability present in every residue, the B-factor, the Eigen value, the variance in both green and red, the elastic network, and the covariance of the complex.

Simulated immune responses indicate an efficient response via antibodies and cytokines
The immunological response that was produced by the in silico vaccine construct was shown to be effective in the simulation results obtained from C-ImmSim.In order to mimic the body's immune response, we gave the vaccine to the animals twice ( Fig 10).The subsequent responses were of a greater quality than the initial ones.The amount of immunoglobulin (IgM) that was created was greater than that of IgG.Following the antigen decrease, the levels of immunoglobulins rose to a high level.This rise could point to the formation of immunological memory as a result of exposure to the antigen, which would show that the individual had practised their immunity.Indicative of these findings was an increase in the quantity of IgM as well as expression of IgG1+IgG2, IgM, and IgG+ IgM.In addition, the upregulation of both cytotoxic T cells and Helper T cells as a result of the multi-epitope vaccine design provided further confirmation of an adequate immune response.
Because of the presence of B-cell isotypes over an extended period of time, isotype switching and the creation of memories in the B cell population are both evidenced.In addition, there was a rise in the number of CTL and HTL cells associated with the formation of memories.In addition to this, macrophage activity increased, although dendritic cell activity remained unchanged.In addition to that, the levels of IFN-beta and TGF beta were both elevated.

Discussion
After an increase in the global health emergence that has been created by various fatal viruses, it is necessary that progress be made in the quest for an acceptable treatment for the diseases.Vaccination is consistently regarded as one of the most successful approaches available for the prevention of infectious diseases caused by viruses [38].The production of vaccinations is widely regarded as one of the most significant steps forward in efforts to cut morbidity and death rates around the world.Not only does it prevent the start of a variety of ailments, but it also identifies a doorway for the disease's removal while simultaneously reducing its toxicity [39].Antigens that are exploited in vaccines do not necessarily need to be virulence factors [40].Despite the fact that most virulence gene products are immunogenic and contribute to acquired immunity protection against disease.Antigens in vaccines can be any infectious disease-causing agent that is known to science.According to Ali et al., protein antigenicity refers to its capacity to activate an immune response to the organism that it is a part of.Antigenicity of a protein is a necessary condition before the protein can be utilized as a component of a vaccine [41].
The present study focused on the identification of the potential epitopes for the development of the effective vaccine using immunoinformatic methods {Bahrami, 2019 64 /id}.In the development of epitope vaccines, the precise identification of epitopes is a crucial stage, and novel approaches to combining epitopes are currently under investigation [42].According to research carried out by Galluzzi et al., the immune response is set off when antigens bind to the ligands or receptors that are located on the surfaces of immune cells [43].It is generally agreed upon that the G and N protein is a necessary antigen for the production of the NiV vaccine.Several vaccines are currently being developed that use NiV G protein or G protein plus F protein as antigens [44][45][46].These vaccines are currently in the process of being tested on humans.It is anticipated that these vaccines will be available in the not-too-distant future.According to the findings, G protein was the primary factor in the production of the antibodies that were capable of neutralizing the virus.The protein subunit platform serves as the basis for the Equivac HeV vaccine, which is the only Hendra virus vaccine that can be used clinically at this time and holds a valid license to do so.This proved that the proposed technology was a trustworthy and established platform, which could be utilized in the production of vaccines against either the Hendra or the NiV [20].
According to a previous study, T-cell epitopes are necessary for MHC molecules to work together and trigger the immune system.When predicting potent T-cell epitopes, it is imperative to select epitopes that are likely to be fixed by MHC [47,48].In addition, determining which cells have CD4+ and CD8+ T receptors is essential for the advancement of MSV [49,50].The incorporation of linkers into specialized sequences has the potential to enhance the construction of vaccines.In the past, many studies have shown that the addition of GPGPG and AAY linkers [41,51,52] between predicted HTL and CTL epitope sequences generated junctional immunogenicity, thereby enabling the logical design and construction of a potent poly-epitope vaccine [53].In their study, Arai et al. [54] utilized the EAAAK linker to connect epitopes and adjuvant, thereby enhancing the bioactivity of the fused protein, achieving a high level of expression, and bolstering the stability of the vaccine construct.Bazan et al. developed a multi-epitope subunit vaccine against the Ebola virus utilizing T-cells.Antigenic epitopes predicted using the Immune Epitope Database (IDEB) were incorporated into the development of a vaccine candidate that exhibited immunogenicity upon expression in rodents [55].Antibodies are able to make a contribution to the eradication of viral loads and the neutralization of pathogen-associated antigens.The ABCpred service was used in order to analyze the antigenic proteome and add epitopes into the design of the vaccine that are capable of activating B cells.Various studies showed the Cholora toxin (CT) adjuvant enhanced the immune response to influenza vaccines in mice and provided cross-protection against heterologous influenza strains and Helicobacter pylori [53,56].Overall, these studies demonstrate the potential of CT as an adjuvant for multi-epitope vaccines against various diseases.CT's ability to induce mucosal immunity and enhance immune responses makes it an attractive candidate for the development of novel vaccines.
The proposed vaccine construct had molecular weight was 63.27 kDa, and its theoretical pI was found to be 9.79, which means it is alkaline.The predicted half-life of our vaccine candidate (30 hours in human reticulocytes, more than 20 hours in yeast, and more than 10 hours in E.) was longer than that of the vaccine made by Waqas et al. [57], which was 7.2 hours in reticulocytes (in vitro) and similar to previous studies [52,58].This means that the vaccine is in contact with the immune system for a long time [59].The vaccine instability index was 28.56, indicating its stability, which was good than some previous studies Kadam et al. 2020; Shafi Mahmud et al., 2021.In general, a protein with an instability index of less than 40 is considered stable [59].The vaccine's aliphatic index was found to be 78.08,which means it stays stable at a lot of different temperatures [60].The aliphatic index value of this vaccine was higher than the aliphatic index of the vaccine designed in works by Sukrit Srivastava, et al  [63].It was found that the multi-epitope vaccine has a GRAVY score of -0.218.This score shows how hydrophilic the vaccine is and is linked to how well proteins dissolve.A negative value means that the vaccine interacts better with water molecules [58].
Understanding the secondary and tertiary structures of proteins is necessary for the production of vaccines, as stated by Meza et al. [53].According to what was discovered in the research, coils make about fifty percent of the secondary structure of the protein.The natural antigenicity of the unfolded native secondary proteins is reflected in the helices, which account for the vast majority of the coils and are composed of helices.The immune system is able to recognize the spiral shape once it has been folded back into its original state.According to the findings of the Ramachandran plot, the 3D model of the vaccine, after being changed, revealed elements that were acceptable and was significantly improved.This is evidenced by the fact that the model was significantly better.The Ramachandran plot reveals that 90.3% of residues begin in preferred domains, whereas only 8.3% of residues begin in allowed regions.This disparity indicates that the output of the model as a whole is satisfactory.The immuno-reactive need in a serological sample, is important for establishing whether or not a vaccine candidate is real.
TLRs are specifically accountable for the detection of viruses and the initiation of innate immunity, which is succeeded by adaptive immunity.Thus, they play a crucial role in innate immunity.TLR3, TLR7, and TLR8 are capable of identifying ds/ssRNA, a prevalent intermediary in the replication process of numerous RNA viruses [64,65].Docking investigations play a pivotal role in comprehending the molecular interactions that regulate immune responses.By simulating the initial recognition stages of the immune system, signaling pathways are stimulated, which ultimately contributes to the development of a robust immune response.In order to enhance the immune response, vaccine adjuvants frequently target TLRs for stimulation.Gaining insight into the manner in which vaccine constructs engage with particular TLRs can enhance the effectiveness of vaccines, customize immune responses, and facilitate the development of personalized vaccines.These studies additionally offer mechanistic understandings regarding the mechanisms by which vaccine constructs stimulate TLRs, thereby improving immunogenicity and informing the development of constructs that stimulate more robust and targeted immune responses [64].In addition, the MD simulation system was used to investigate the connection of TLR on immune cells via the vaccine-TLR docking complex.This investigation was carried out in order togather more information.According to the results of our study, TLRs displayed a strong affinity for the newly designed vaccine.According to Crozat and Beutler [66], the interaction between TLR and the vaccination revealed that it had the potential to generate both innate and adaptive inflammatory responses.LigPlot+ was used to illustrate the intramolecular interactions between the vaccine design and TLR3, 7 and 8.The effectiveness of the vaccine design to induce an immune response is explained by an immunological simulation study.One of the requirements for becoming an effective vaccination candidate is the induction of memory B-cells and T-cells [67,68].The number of memory B-cells and T-cells grew with each administration of the vaccine construct, and this population level persisted after the third injection.The ability of vaccine designs to create an antiviral state is also demonstrated by the raised levels of IFN-γ and IL2.In general, IFN-γ participates in antiviral replication and serves as the primary effector molecule in cell-mediated immunity.Escherichia coli is the best host to use if you want to produce a larger amount of recombinant vaccine protein.Java The positive content of vaccine constructions for high-level protein production in the E. coli host was confirmed by the codon adaptability tool, which showed that the codon adaptability index for vaccine candidates was > 0.95% and that the GC content was 65%.
The traditional method of developing vaccines is an expensive and drawn-out process that necessitates growing the pathogen in a lab [69].Without cultivating the pathogen, reverse vaccinology uses computational analysis and genomic data to generate vaccines.Predicting T and B cell epitopes, antigen processing, antigenicity, allergenicity, toxicity prediction, and TLRpeptide docking are crucial steps in designing virus vaccines [70].

Conclusion
Immune informatics tools were used in this study to develop a multi-epitope NiV vaccine.There have only been a few studies that have identified epitopes that target T-cells and B-cells using NiV's attachment glycoproteins (G) and nucleocapids (N).The present NiV-G and NiV-NP, MHC I and MHC II epitope vaccine formulations are highly immunogenic and antigenic.The adjuvant was added to further increase the construct's efficacy after the epitopes were coupled using the appropriate linkers.The ability of the chimeric vaccine peptide to elicit an immune response must be tested in vitro and in vivo.It would facilitate the future creation of NiV vaccine that would benefit the intended population.

3. 7 .
Prediction, validation and refinement of the tertiary structure I-TASSER predicted a total of 5 models utilizing 10 threading templates with Z-scores, and the models were chosen due to their high confidence C values, which vary from -0.63 to -4.18.Five models were created using the GalaxyRefine server to refine the chimeric vaccine's 3D structure (Fig 5A).Model 1 was chosen for molecular dynamic (MD) simulation based on its model quality score, which included GDT-HA (0.9318) and RMSD (0.480) (Fig 5B).The 3D structure was validated to look for quality issues and probable mistakes.The 3D structure analysis performed by ProSA-web after the MD simulation shows that the structure had a Z-score of 7.52 or less (Fig 5C).According to the Ramachandran plot generated by PROCHECK, 90.3% of the residues are in the regions with the most favorable conditions, 8.3% are in the further permitted regions, and 0.2% are in the liberally permitted regions with no residue in the prohibited regions (Fig 5D).

Fig 5 .
Fig 5. Structure prediction, refinement and validation of the 3D structure of the prepared vaccine construct (a) predicted 3d structure by I-TASSER (b) refinement of the 3D structure of the chimeric vaccine by using galaxyRefine server (c) ProSA Z-score plot shows a -7.52 score in the range of conformation of the predicted protein (d) PROCHECK's Ramachandran plot illustrates that the residues are placed in the allowed and favored regions.https://doi.org/10.1371/journal.pone.0300507.g005 Fig 6A), TLR 7 (Fig 6B) and TLR 8 (Fig 6C) and vaccine construct were analyzed by LigPlot software.

Fig 6 .
Fig 6.The docked complex of the multi-epitope vaccine and TLR was visualized using the UCSF Chimera software and interaction in the docked complex has been shown using Ligplot software.(a) docked complex and interaction for TLR3 with vaccine construct (b) docked complex and interaction for TLR7 with vaccine construct (c) docked complex and interaction for TLR8 with vaccine construct.https://doi.org/10.1371/journal.pone.0300507.g006

Table 2 . MHC-I epitope prediction for the nucleocapsid protein.
Population coverage of the selected epitopes was 84.51% in Central Africa, 90.22% in East Africa, 89.16% in North Africa, 93.63% in South Africa, and 94.25% in West Africa, respectively.The percentages for East Asia, Northeast Asia, South Asia, Southeast Asia, and Southwest Asia were 96.54%, 92.1%, 89.45%, 92.1%, and 88.6%, respectively.Fig 2 shows that the iedb.org/population/.Fig 2 shows that the MHC I and MHC II epitopes that are employed for population coverage study may cover 96.76% of the world's population.

Table 5 . List of linear B-cell epitopes with their location and score by using the ABCpred server in NIPAH glycoprotein. Rank Sequence Start position Score
The final 572 amino acid peptide vaccine design (Fig3) was created.The physiochemical parameters of the generated 572 amino acid long vaccine construct were determined using the protparam server, which revealed a molecular weight of 63274.29 and a theoretical pI of 9.79, indicating its basic nature and computed instability index of 28.56 indicates that the construct is stable.The aliphatic index was 78.08, which indicates that the substance is thermostable throughout a range of temperatures.Additionally, the grand average of hydropathicity (GRAVY) was -0.218, which shows that the vaccine design is hydrophobic.Protein-Sol server score of 0.608 indicated that our vaccine construct had a good solubility rate.

Table 6 . List of linear B-cell epitopes with their location and score by using the ABCpred server in NIPAH Nucleocapsid protein.
https://doi.org/10.1371/journal.pone.0300507.t006